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We study the properties of a non-Gaussian density matrix for an O(N) scalar field in the context 
of the incomplete description picture. This is of relevance for studies of decoherence and entropy 
production in quantum field theory. In particular, we study how the inclusion of the simplest non- 
Gaussian correlator in the set of measured observables modifies the effective (Gaussian) description 
one can infer from the knowledge of the two-point functions only. We compute exactly the matrix 
elements of the non-Gaussian density matrix at leading order in a 1/iV-expansion. This allows us to 
study the quantum properties (purity, entropy, coherence) of the corresponding state for arbitrarily 
strong nongaussianity. We find that if the Gaussian and the non-Gaussian observers essentially 
agree concerning quantum purity or correlation entropy, their conclusion can significantly differ for 
other, more detailed aspects such as the degree of quantum coherence of the system. 
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I. INTRODUCTION 

To completely specify the state of a (quantum) system 
requires one to perform as many measurements as there 
are independent degrees of freedom (d.o.f.). This may be 
difficult in practice for systems with a large number of 
d.o.f. or even impossible when the latter is infinite, e.g. 
for (quantum) fields. In fact, one often has to infer the 
effective state of the system from a restricted amount of 
information, which may or may not give a good descrip- 
tion of the actual state of the system. 

This incomplete description picture [l[ provides a ba- 
sis for various studies of entropy production in quantum 
field theory (QFT) 0,[|. In this context, one defines an 
entropy which measures the amount of missing informa- 
tion in the set of measured observables. More recently, 
similar ideas have been advocated to discuss quantum 
(de)coherence in QFT 0-16}. Both issues are of great 
interest in various contexts such as inflationnary cosmol- 
ogy ! baryogenesis [1] , neutrino [1] and quark-gluon 
plasma physics |l0| , or the physics of ultracold quantum 
gases jllj . 

First principle calculations of the real-time dynamics 
of entropy production and quantum decoherence in QFT 
have been performed recently in Refs. [5|, lfj( in the case 
where the set of measured observables is restricted to 
Gaussian, two-point, correlators of the fields. As a result 
of the nonlinear evolution, information flows from the 
subset of Gaussian correlations to - unmeasured - higher 
correlations. In it was demonstrated in the context 



of a self-interacting O(N) scalar field that, even in the 
absence of an environment, a pure, coherent quantum 
state eventually appears as mixed and decohered, and 
at late times even as a thermal state, to the observer 
restricted to Gaussian correlators - we shall call the later 
the Gaussian observer - even though the time evolution 
is unitary. The case of a field coupled to an environment 
has been considered in [(|, with similar conclusions. 

Very interesting questions in this context concern the 
inclusion of higher correlation functions in the set of mea- 
sured observables. How much information one recovers 
by measuring the first non-Gaussian correlators? Is the 
spread of information in the space of correlation func- 
tions uniform? Or is it mainly concentrated in a subset 
of relevant observables? Besides the above-mentionned 
applications, such questions are of general interest for 
the understanding of nonequilibrium dynamics and ther- 
malization in QFT [Hj]. 
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In [13j Koksma, Prokopec and Schmidt have studied 
in detail the question of the correlation entropy for var- 
ious deformations of Gaussian density matrices in scalar 
field theories. Here, we consider the non- Gaussian den- 
sity matrix one obtains in the incomplete description pic- 
ture by including the simplest non-Gaussian correlator - 
namely the field four-point function - in the set of mea- 
sured observables. We study how the inclusion of this bit 
of information modifies the conclusions one would infer 
from the knowledge of the two-point functions only. 

We find that, depending on what they ask, the Gaus- 
sian and non-Gaussian observers may arrive at very dif- 
ferent conclusions. In particular, if they essentially agree 
concerning global properties such as the quantum pu- 
rity or the entropy of the system, whatever the strength 
of the nongaussianity, they can be led to very different 
conclusions concerning other, more detailed aspects of 
the quantum state, such as the degree of quantum co- 
herence or the shape of the probability distributions of 
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various observables. For instance, we exhibit situations 
where the Gaussian observer concludes to an essentially 
thermal state, whereas the non-Gaussian observer sees a 
highly coherent quantum state. 

Specifically, we consider a situation with space trans- 
lation invariance, where the fields can be described by 
their Fourier modes and we assume our observers mea- 
sure Gaussian and non-Gaussian correlators of separate 
modes. In this case, the effective density operator is a 
direct product of individual density operators for each 
Fourier mode. We focus here on one of these individual 
operators - that of the zero mode to be precise. This 
is essentially equivalent to considering the density oper- 
ator of a quantum mechanical model with one degree of 
freedom. 

In principle this could be computed exactly with nu- 
merical techniques. Instead, we want to have as much 
analytical understanding as possible as well as simple 
practical tools which can be easily implemented in fu- 
ture nonequilibrium calculations [14]. For this purpose, 
we consider an TV-uplet of scalar fields in an O(N)- 
symmetric state and compute the desired properties in 
the limit N — ¥ oo, for which we get semi-analytical re- 
sults. Moreover, such a nonperturbative approach allows 
us to study the case of arbitrarily large nongaussianities. 

We describe the reduced density matrix inferred by 
the "non-Gaussian" observer, in Section |TTJ Some global, 
basis-independent, properties of the corresponding quan- 
tum state, such as correlation entropy or quantum purity, 
are analyzed in Section IIIII In Section IIV1 we present 
the calculation (new to our knowledge) of the matrix el- 
ements of the non-Gaussian density operator in the limit 
N — » oo. We analyze the detailed properties of the quan- 
tum state, compute the Wigner function on phase space, 
and study the issue of quantum coherence in Section fVl 
A number of appendices are devoted to technical details. 



II. THE REDUCED DENSITY MATRIX 

We consider N pairs of canonically conjugate scalar 
field operators (ip ai % a ),a = l,...,N. We assume a 
0(-/V)-symmetric state with given Gaussian correlators 

(TTaTTb) = S ab K, (tp a (pb) = S ab F, (lfi a TT b + TT b (f> a ) = 2S ab R, 

(1) 

and non-Gaussian four-field correlator 



(ip a ip b ip c lfi d ) = (6 ab S c d + Sac5 b d + 5adSbc) 



F + 



N 



(2) 



Here C4/N denotes the non-trivial, connected contribu- 
tion. The least biased effective density matrix D one can 
infer from the knowledge of these correlators is the one 
which is consistent with the measured observables, i.e. 
(f 2 ) = tr[D ip 2 } etc., and which maximizes the amount 
of missing information, measured by the von Neumann 
entropy S = — tr[Z?ln£)]. It is given by [l[ 

D = 4 cxp(-J") with Z = trexp(-J') , (3) 



where (we extract an explicit factor TV for the purpose of 
the 1 /TV-expansion) 

T = An a n a +Btp a tp a +C(tp a n a +7r a ip a ) + ^((p a fi a ) 2 . (4) 

The parameters A, B, C and 77 must be adjusted so that 
the correlators K, F, R and C4 obtained from D agree 
with their measured values. For instance, they can be 
obtained from the following relations 

d A lnZ d B hxZ d c lnZ ... 
K = — — , t 1 = — — , K = — (5) 



N 



and 



N 



F 2 , _ d n In Z 



2N 



N N + 2 



A. Large-TV limit 



(6) 



Following standard procedures, the partition function 
Z can be given the path integral representation 

Z = I X> ( p2? 7re /o ldT {^a(-r)7r a (T)-J r ( V (r),ir(T))} fj^ 

J periodic 

where up = dip/dr and where the functional integral is to 
be performed on periodic configurations f a {0) = </?a(l) 
[2^ ]. Performing the n integration, one gets 

Z = Af(A N ) [ Vtpe^i*"*^- 3 ')^-*^ 2 } 

J per. 

(8) 

where we used the fact that, for periodic field configu- 
rations, fdnpatpa — 0. Here B' = B — C 2 j A and, up 
to an infinite A-independent normalization factor (note 
that Af(AB) = M(A)M(B)) 

AT (A) = fv X e- dT A ^ oc e" ^ Tr ln A (9) 



where Tr denotes a functional trace. 
Making use of the standard trick pi 



Vxe 



(10) 



we can perform the (^-integration to get, up to an irrele- 
vant factor 

where 

G-\x; r, r') = (~ + AA[B' i X {r)]\ S(r r'). 

(12) 

This is the starting point for the standard 1 /^-expansion 
[l6l |. The limit N — > 00 is given by the saddle point 
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approximation: One has, up to an irrelevant additive 
constant, 



— = --TrLnG ( X ) - - , 

where \ is given by the saddle-point equation 
8Tr LnG~ 1 (x) 



(13) 



n- 



Sx(r) 



AiAr)TrG(x). (14) 



For constant x, the functional trace on RHS can be evalu- 
ated by standard techniques (see the details in Appendix 
IB"]) and Eq. (|14|) can be rewritten as 



1 



ztanh(z/2) ' 



(15) 



where £ = 8A 2 r], z Q = 2VAB', z = 2y/A(B' - igj. We 
mention that Zq plays the role of a mass squared-like pa- 
rameter in the Gaussian case, governing the width of the 
(/^-distribution and is constrained to be positive. In the 
non-Gaussian z 2 plays the role of a renormalized mass 
parameter. In Appendix [Cl we show that only real solu- 
tions, i.e. z 2 > 0, are physically allowed. One easily sees 
that Eq. (|15j) always admits one and only one real solu- 
tion, which can be obtained numerically. Note that one 
recovers the Gaussian result z 2 = z 2 in the limit r\ — > 0. 
Notice also that, for r} ^ 0, z 2 can be either positive or 
negative. 

Expressions of the correlators K, F, R in the limit 
iV^co can be directly obtained from (fl"3)) . using ([5]). 
The connected correlator C4/N is 0(1/N) and necessi- 
tates a more involved calculation. The latter is detailed 
in Appendix [A"l In the limit N — > 00, one obtains the 
parameters of ((4]) in terms of the correlators as [13] 



A = kF , B = kK — 2r]F , C 



-kR . 



(16) 



where 



ln(l + l/n) 
2n + l 



and n + - = \J FK - R 2 . (17) 



It is remarkable that only the parameter B gets mod- 
ified by nongaussiantiy. This is a consequence of both 
our choice of non-Gaussian correlator and the large-TV 
approximation. Note that n = (e z — 1) , with z the 
solution of the gap equation above. The value of r\ can 
be obtained from the relation (see Appendix IB"]) 



1 



ln(l+l/n) 



c 2 



1 + C,x 1 + x 



where we defined 



V ( F 
C = 1 + 2Kn(n + 1) and x = — 



k \n + 1/2 



(18) 



(19) 
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FIG. 1: The ratio C 4 /2F 2 as function of x for n = (blue), 
n — 1 (green) and n = 10 (red). The red and blue curves 
are well described by the limiting functions (|22[) and (|23[) 
respectively. 



and (notice that /(1/2) 
1 



2n + l) 



m = 



1 (l+nyy + n 



,2y 



2ytanh[j/ln(l + l/n)] 2y (l + n) 2 « - n 2 ^ 



(20) 

In the perturbative - near Gaussian - regime, small val- 
ues of the connected correlator C4 as compared to the 
disconnected contribution F 2 corresponds to small val- 
ues of x. In this regime, one gets the linear relation 



g 4 
2F 2 



2{2n 



(l + 2n(n + l)(3C + 2)). (21) 



In contrast, for arbitrarily large values of x, the ratio 
C4/F 2 saturates due to non-perturbative effects. Wc 
see that —1 < C4/2F 2 < 0. It is interesting to no- 
tice that the extremal value C4 = — 2F 2 actually ex- 
actly cancels the corresponding (in the l/A^-expansion) 
connected contribution to the correlator ((ip a <p a ) 2 /N) — 
NF 2 + (2F 2 + C 4 ) + 0(1/N). To separate between 
the perturbative (near Gaussian) and non-perturbative 
(strongly non-Gausssian) regimes, we note that the Gaus- 
sian square-mass parameter Zq = ln 2 (l + l/n)(l — 2x) 
turns negative for x > 1/2. We expect a qualitative 
change in the properties of the density operator when 
this happens. 

We note that, as a function of x the ratio C4/2F 2 is 
bounded from below by the large-n limit: 



C 4 
2F 2 



2x 



n»l 



1 + 2x 



and from above by the small-n limit: 

C 4 
2F 2 



n<l 



1 



(22) 



(23) 



We plot the ratio C4/2F 2 as a function of x for various 
values of n in Fig. Q] 

Finally we mention that, using x = 2irjF ', one has 

lnZ 



In y/n(n + 1) + -rK (2n + 1) (24) 
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The first term is the standard Gaussian contribution and 
the second one the non-Gaussian correction. We now 
turn to the analysis of the properties of the density ma- 
trix ©. 

III. CORRELATION ENTROPY, QUANTUM 
PURITY 

For the Gaussian observer n is the only intrinsic - ba- 
sis independent - parameter which characterizes the state 
of the system. Any basis independent observable, such 
as correlation entropy or quantum purity, is fully deter- 
mined by the value of n. In this section, we analyze this 
kind of global intrinsic observables in the non-Gaussian 
case. 
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FIG. 2: The ratio of non-Gaussian vs. Gaussian purities p/po 
as a function of the nongaussianity x for, from top to bottom, 
n = 0, 0.1, 0.5, 1 and 10. The latter case is well-described by 
the n > 1 limit, Eq. ((33)) . 



A. Correlation entropy 

The correlation entropy 0, 0] measures the amount of 
missing information, i.e. the amount of information not 
contained in the subset of measured observables. It is 
defined as 



S = -tr(D In D) = In Z + (T) 



(25) 



and is easily evaluated in the large- N approximation, us- 
ing (g]), (pi]) and l[2"3|). One finds that the contribution 
from non-Gaussian terms exactly cancel and one is left 
with the Gaussian expression [28[ : 



- = (n + 1) ln(n + 1) 



i In i 



(26) 



This shows that the inclusion of the correlator C4 does 
not bring any information, in the sense of that measured 
by entropy. This may seem surprising at first sight since 
there is clearly some "information" in the nongaussianity. 
But in fact, entropy is a one number which measures 
the global amount of (missing) information, . There can 
clearly be an infinity of different density operators having 
the same entropy. 



B. Quantum purity 



The quantity P = tvD 2 < 1 measures the quantum 
purity of the state described by the density operator D. 
It can be computed as 



P = trLH = 



tr e 



-2T 



(tr e 



-^2 



Z{2A,2B,2C, 2rj) 
Z 2 (A,B,C,r 1 ) 



(27) 



where the numerator is obtained as in the previous sec- 
tion by replacing all parameters by twice their values. 
Writing the solution of the corresponding saddle-point 
equation z(2A, 2B, 2C, 2r\) = 25, the latter can be writ- 
ten 



; tanh ; 



(28) 



In the large- N limit, the N field components are essen- 
tially identical and independent degrees of freedom and 
the total purity can be written as P — p N where p can be 
viewed as the individual purity of each d.o.f. [29j. Intro- 
ducing h — (e z — 1) _1 and using Eq. p4|) . we get, after 
some calculations, 



1 



P 



fl(h + 1) -I KK (2n+l) 2 {] 



2h + 1 n(n + 1) 



LT+4ST 



(29) 

where k = ln(l + l/n)/(2h + 1) = 2tanh(5/2). Notice 
that, for real positive x, tanh a; > tanh(x/2) and, there- 
fore, z < z, which implies n > n and k > k. Notice also 
that purity only depends on the parameters n and x. 

We check that we recover the usual result [5J in the 
Gaussian case (x = 0), for which z — z = zq: 



Po 



2n+ 1 



(30) 



The pure Gaussian state corresponds to n -C 1 . A simple 
analysis shows that the case n^Cl^z^lis equivalent 
to n < 1 i > 1. In this situation, one gets, from Eqs. 
([151) and (|2"g]). z w z and thus 



p w p rj 1 for n -C 1. 



(31) 



We conclude that both the Gaussian and the non- 
Gaussian observers agree on the purity of the system 
when the latter is pure. 

The other extreme is n > 1, which implies ti>l. In 
this case, one has zwl/nCl and z rj 1/h 1. The 
saddle point equation (|28[) simplifies to a quadratic equa- 
tion for z 2 and one gets, after some simple calculations: 



~ 2 
n 



l-2x+ VI + 4x 2 



and 



P_ 

Po 



-xfl-- 
2 V 



for «>1. 



(32) 



(33) 
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One easily checks that 1 > p/po > \/2/e ~ 0.86, the 
two limits corresponding to the Gaussian (x = 0) and 
the extreme non-Gaussian (x 3> 1) cases respectively. 
We conclude that both observers essentially agree on the 
degree of quantum purity of the system whatever the 
value of n and the strength of the nongaussianity. This 
is demonstrated on Fig. [2] which shows the individual 
purity p as a function of the nongaussianity x for given 
Gaussian correlators, i.e. for fixed n. 

Thus we see that both the correlation entropy and the 
quantum purity are pretty insensitive to the degree of 
nongaussiantiy of the system [3(J. Such global observ- 
ables are not sensitive to high moments of the under- 
lying distribution and thus do not contain enough infor- 
mation to qualitatively distinguish between the Gaussian 
and non-Gaussian cases. 

Still, as explained above, we expect substantial, qual- 
itative changes in the detailed properties of the density 
operators inferred by our two observers. In the next sec- 
tion, we compute the detailed structure of the density 
operator in the large- N limit and study how it evolves as 
a function of the nongaussianity. 



IV. MATRIX ELEMENTS IN THE LIMIT N -> oo 

We now wish to compute the matrix elements of the 
density operator in various basis. The "position" basis, 
i.e. the basis of eigenvalues of the field operators, turns 
out to be particularly suited to take the N — > oo limit. 
We start with the general functional integral formula: 



(<P2\e -H^i) 



T>(pT>n e 



Jo dT {iVaKa-Fip,*)} 



(34) 



where \<p) = J\ \<Pa)- Here the functional integral is to 
be performed on field configurations such that <p(0) = tpi 
and <p(l) = 2, where arrows denote O(N) vectors. Note 
that, because of O(N) symmetry, the matrix element (|34p 
only depend on the invariants ip\ , ip\ and 0i -02- It is 
useful to introduce the rescaled quantities [31 1 



NA 



ANA 



and w 



if ■ s 
2NA 



where 



f = 



2 + Cp x 



and s = tfi2 — <fii 



(35) 



(36) 



Using Eq. ([ID]), we write 



(ip 2 \e ^Itpi) = TV 



where we defined 



V X e~ N ( S * 2 ' n[x]+S ° dT ^) , 
(37) 



Jtpx 



(A,B-ix,C)} 

1 

(38) 



with the quadratic form 

7"q(A-B)C) = A-K a Tr a +Btp a (p a + C(tp a -K a + % a tp a ). (39) 

As for the calculation of the partition function above, the 
large- iV limit is given by the saddle-point contribution: 



(He ^IVl) 



^ q (A,s- l x(t l 2 ,t, 2 ),c)| V5 ^ e -i|x 2 (« 2 ^ 2 ; 



= (He 

where the saddle-point x(u 2 ,v 2 ), given by 



\(u 2 ,v 2 ) SS Va , Vl \x] 



2>] 



,(40) 



(41) 



x(t)=x(« 2 ^ 2 ) 



only depends u 2 and v 2 defined in Eq. (|35j> . as shown 
below. 

We note that, using (l38|) . the saddle-point equation 
(|4ip can be rewritten as 

(42) 

Therefore, we find that, in the limit N — > oo, the ma- 
trix elements of the non-Gaussian operator exp(— J 7 ) can 
be entirely expressed in terms of those of the Gaussian 
operator exp(— J- q ). The latter are well-known: 

(He~^ q |</3i) = J^e _A, { Fo(zo)+Fu(zo)ll2+F " (zo)l ' 2+2lC,i, J', 
A% 

where zq has ben defined in Eq. (|15[) and u,v,w in 
Eq. ([33)) and where 



Fo(z) 
F u {z) 
F v {z) 



In 



4-7T sinh(z) ' 



Z ( z 

— tanh - 
2 V2 



z/2 



tanh(z/2) 



(44) 
(45) 
(46) 



We emphasize that Eqs. (gOJ and (|42J) in fact hold 
in any basis and, more generally for any linear combi- 
nation of matrix elements (32j. For instance, Eqs. (|13[) 
and (TT41 can be obtained as a particular case when one 
replace (<^2 1 e — | </?i ) — > tre^^i in both Eqs. (|40|) and 
(|42p . Similarly, we mention that similar equations can 
be obtained for the Wigner function ([6"0"]) to be discussed 
below. 

The position basis employed here presents the advan- 
tage that the solutions of Eq. (|42p are particularly sim- 
ple. Notice first that the phase of (<f2\e~' Fq \«Pi) in (O 
is the only to-dependent term and is independent of B. 
Thus this term does not contribute to the saddle-point 
equation, which implies, first, that zx(u 2 ,i> 2 ) is real and, 
second, that it only depends on u 2 and v 2 . In terms of 
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U 



FIG. 3: The function d(u 2 ,v 2 ), normalized to its maximum 
value in the (u, n)-plane for n = 10 and x = 15. 



U 

FIG. 4: The function d(u 2 , 0) as a function of u for fixed Gaus- 
sian correlators with n — 10 and increasing nongaussiantiy x: 
from the Gaussian case x = (blue) to x — 0.5 (green), x = 1 
(orange) and x = 15 (red). The corresponding values of the 
ratio \C 4 /2F 2 \ can be obtained from Eq. J22J: 0, 0.5, 0.67 and 
0.96 respectively. The red curve is the v = slice of Fig. O 



the variable z(u 2 ,v 2 ) = A[B' — ix(u 2 , v 2 )] we obtain 
the following equation 

2 2 

- = fo(z) + fu(z)u 2 + f v (z)v 2 (47) 



where fi(z) = 2F-(z)/z, with i = Q,u,v; Explicitly: 



fu(z) = 



fv(z) = 



z 2 \tanh(z) 
1 

2cosh 2 (z/2) 
1 



(48) 



sinh(z) \ 

— + lj , (49) 

sinh(z) 



- 1 



(50) 



2sinh 2 (z/2) V z 
In terms of the solution z = z(u 2 , v 2 ) of Eq. (|4T1) . we have 

( W! |i?|v»i>=A-T2J( U 2 ,« 2 )e- aw<7 "' > (51) 
where I?(u 2 ,w 2 ) is a real function given by 

(±41 



V(u 2 ,v 2 ) = -e I 

(52) 

We seek solutions of Eq. (jT7|) with z 2 real, i.e. z can be 
either real or purely imaginary. We show in Appendix [Cl 
that physically allowed solutions are such that z 2 > — it 2 . 
It is rather simple to see, e.g. graphically, that Eq. (|T7)) 
always admits one and only one such solution. Moreover, 
we readily see from the small z behavior of the RHS of 
Eq. 63), namely / (0) = /„(0) = 1/3 and /„(0) = 1, 
that Eq. (jTT)) admits a real solution (z 2 > 0) if 



IT . Zn 



- + U -\ > 

3 3 - e 



(53) 



and a purely imaginary one (— 7r 2 < z 2 < 0) otherwise. 
A condition for the existence of imaginary solutions is, 
therefore, Zq/£ < —1/3. Noticing that 



1/1 



k \2x 



1 



this condition becomes 



K < 3 and x > 



1 



2(1 - k/3) 



(54) 



(55) 



In particular, we see that no imaginary solution can ap- 
pear, whatever the strength of the nongaussiantity x if 



k > 3, which is equivalent to n < n c 



We also 



see that, for n > n c , the appearance of imaginary solu- 
tions is only possible for large enough nongaussianity x. 
In particular, this necessitates that z 2 /^ be sufficiently 
negative and is, therefore, a nonperturbative effect. 

When conditions (|55|) are satisfied, there is a region of 
the (w, v) plane - for sufficiently small u 2 and v 2 - where 
z = z(u 2 ,v 2 ) is purely imaginary. Writing 



v(u 2 y) = d N {u*y) 



and noticing that 



and 



d u 2ln d(u 2 ,v 2 ) = -F u (z). 



d v 2\nd(u 2 , v 2 ) = —F v (z) 



(56) 



(57) 



(58) 



and that, for z 2 > —tt 2 , sign(F u (z)) = sign(z 2 ) and 
F v (z) > 0, we conclude - see Appendix [E] - that 
lnT>(u 2 ,v 2 ) is always a monotonously decreasing func- 
tion of v whereas it either is monotonously decreasing in 
the it-direction if condition (|53[) is not satisfied or if v 
is large enough: v 2 < — z 2 /£ — 1/3, or has a maximum 
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at u = u c (v 2 ) otherwise, corresponding to the condition 
z(u%v 2 ) = 0, i.e. 



(59) 



This is illustrated on Fig. [31 where we plot the function 
d(u 2 , v 2 ) as a function of u and v for Gaussian parameter 
n = 10 and non- Gaussian one x = 15. We see that when- 
ever condition (1551) is fulfilled, the non-Gaussian distri- 
bution qualitatively differs from the corresponding Gaus- 
sian one. We show, in Fig.[4]the v — slice of the function 
d, that is essentially the probability distribution of field 
values ((p\D\ip), for fixed Gaussian parameter n — 10 and 
various nongaussianties x, ranging from x = (Gaussian) 
to x = 15 [33]. We mention that similar shapes have been 
also obtained in Ref. [l3j . 

As expected, the structure of the density operators in 
field-configuration space inferred by either the Gaussian 
or the non-Gaussian observers significantly differ for large 
enough nongaussiantity. In the next section we further 
illustrate this difference by computing the Wigner distri- 
bution on phase-space and studying the quantum coher- 
ence properties inferred by both observers. 



V. WIGNER FUNCTION AND QUANTUM 
COHERENCE 



For a Gaussian density operator the Wigner function 
(see below) is positive definite and can be interpreted as 
a probability distribution on phase-space. This remains 
true for the non-Gaussian density operator Q in the 
large- N limit. One interesting property one can directly 
read on the Wigner function is the degree of quantum 
coherence of the system in the coherent-state basis. In- 
deed, if the Heisenberg uncertainty principle guarantees 
that the area of phase space where the Wigner function is 
non-vanishing cannot be smaller than 1 (in appropriate 
units), one can stil have a very squeezed distribution in a 
given direction. As recalled in Appendix [Dl this gives rise 
to non-trivial quantum coherence between distant semi- 
classical (coherent) states [34J. We analyze the shape of 
the Wigner function for our non-Gaussian density oper- 
ator in this section. 



A. Wigner function 

The Wigner function is defined as 
W(0, 7?) = J d N s e ijf - s (<p + s/2\ D\ip- a/2) , (60) 
where 

(<p + s/2\D\<p-s/2) = A-^V{u 2 ,v 2 )e- l ^\ (61) 

with T>(u 2 ,v 2 ) real. Eq. (|61|) and O(N) symmetry imply 
that the Wigner function really depends on two- variables 




FIG. 5: The wigner function w(u 2 , r 2 ), normalized to its max- 
imum, for the same values of parameters as in Fig. [3] 



only: 



W(ip,Tf) = W(w 2 ,r 2 ), 



where we defined 



AA 

~w 



'„ c _ 



(62) 



(63) 



Exploiting the O(N) symmetry, N — 2 angular integra- 
tions can be trivially performed in Eq. (1601) and the re- 
maining angular integral is given by a Bessel function 
J v (x). One obtains 



W{u 2 



Nr 



dvv^ jN_ 1 (Nrv)V(u 2 



(64) 

It is easy to check that, at large N, the integral is dom- 
inated by values v ~ 1 due to the combination of the 
phase-space factor v N / 2 and of the rapid decrease of the 
function T> at large v. It follows that for the typical values 
of r we are interested in, that is r ~ 1, the argument of 
the Bessel function is of the same order as its index ~ N, 
which is precisely the region where the Bessel function 
cannot be given a simple approximation. We com put e 
the large- N limit of the integral numerically |35| . 
Fig.[5]shows the function w(u 2 ,r 2 ) defined as 



W{u 2 y)=w N {u 2 ,v 2 ) 



(65) 



as a function of u and r for the same values of parameters 
as in Fig. [3J 

As for the function d(u 2 ,v 2 ) before, the shape of 
w(u 2 ,r 2 ) only depends on n and x. This encodes the 
Wigner distribution on the A f2 -dimensional phase-space, 
the translation to which depends on the actual values of 
the Gaussian correlators F and R through Eqs. (HI 



(|35| and ((63)) . Here we illustrate the content of the func- 
tion w(u 2 ,r 2 ) in the cases where the O(N) vectors <p and 
7? are either collinear or perpendicular to each other. It 
is useful to use the parametrization of Gaussian correla- 
tors introduced in [5j (up to a redefinition of the angle 
4> — > vr — 4>): 

F = a (1 + 7Cos<^>), 

K = a (1 — 7cos0), (66) 

R = a 7 sin (p, 

with a = (n + 1/2)/ y/l - 7 2 , < 7 < 1 and < cj> < 2tt. 
The parameter 7 controls the overall squeezing of the 
Wigner distribution. Note that C/A = R/F = tan(0/2). 

When ip and 7? are collinear, one has r 2 oc (w — 
95 tan(</>/2)) 2 . The angle gives a positive y? 2 term, pro- 
portional to F, which adds to the already present one 
and thus squeezes the distribution in the ^-direction. It 
also contributes a (pit-term which rotates the main axis 
of the distribution by an angle </>/2, as can be seen on 
Fig. [6] In contrast, in the case where (p and if are per- 
pendicular, one has r 2 cx n 2 + <p 2 tan 2 (<rj/2), which results 
in a squeezing of the distribution in the (^-direction (by 
the same amount as before), but there is no ipn-term and 
thus no rotation. The case = corresponds to a over- 
all squeezing (for 7 ^ 0) in the 7r-direction and an overal 
stretching in the (^-direction, i.e. F > K, whereas the 
case <j> = 7r gives the opposite effect. 

We already see from Fig. [5] that the Wigner distri- 
bution can be strongly squeezed in some direction and 
therefore present a high degree of quantum coherence. 
We discuss this further in the next section and compare 
to the Gaussian case. 



B. Quantum coherence 

To illustrate the implications of the very different 
Gaussian vs. non-Gaussian Wigner functions for quan- 
tum coherence, we consider the case where the correlator 
R = 0, that is r 2 — AA-k 2 /N . As described before, for suf- 
ficiently large values of n and x [36[ , the function d(u 2 , v 2 ) 
exhibits a clear peak-structure in the (u, u)-plane, located 
at u — Uq, v = 0, with, see Eq. (|59|) . Uq = —z 2 j^ — 1/3. 
We find that for x ^> n 2 >• 1 the peak is well-described 
by a Gaussian: 

d(u 2 ,v 2 ) « d(ul,0)e (67) 
with, sec Appendix IB"! 

1 n 2 1 
u 2 = 2n 2 , 8l = - + — and S 2 V = -. (68) 

In that case, the Wigner function is also Gaussian: 

w(u 2 ,r 2 ) w w(ul,0)e 2S 'i (69) 
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-10 -5 5 10 -10 -5 5 10 
q)/V /2 

FIG. 6: Contour plots of the Wigner distribution w in physical 
coordinates for the cases where (p and 7? are either parallel 
[para) or perpendicular (perp) in O(N) space, for the same 
parameters as in Fig. [5] We use the parametrization (|66|) 
with n — 10 and 7 = 0.9. It should be kept in mind that the 
above figures are two-dimensional slices of a 0(iV)-symnietric 
distribution in iV 2 -dimensional phase-space. In particular, 
the two apparently separated structures in each figure are in 
fact connected by 0(iV)-transformations. 

with 

S 2 . = 2. (70) 

In terms of physical variables ip and 7T, the distribution 
peaks at (p^/N = Au^, with widths S 2 /N = AS 2 and 
51/ N = S 2 /4:A, where A = kF « F/2n 2 w 1/2K, that 
is: 

4-f, S = _^ and & = K . 

N ' N Yin 2 N v ' 

Note that the width in the 7r-direction is unchanged com- 
pared to the corresponding Gaussian (x — 0) case. In 
contrast, the width in the (^-direction is much smaller 
than the corresponding Gaussian one ~ F [371 ] . 

We see that our observers may arrive at very different 
conclusions concerning the degree of quantum coherence 
of the system. Indeed, suppose that F ~ K ~ n, in which 
case the Gaussian observer concludes to a statistical mix- 
ture of many uncorrelated semi-classical states, see e.g. 
[3,0]. In constrast, the non-Gaussian observer gets an 
0(7V)-symmetric, highly squeezed Wigner distribution 
around a non-zero value tp — ipo, with S V /\^N ~ 1/n « 1 
and thus concludes to a high degree of quantum coher- 
ence (see Appendix |P|l . This situation is illustrated in 

Fig.ra 

Also interesting is the case where F ~ 1 and K ~ 
n 2 , for which the Gaussian observer concludes to what 
Campo and Parentani call a partially decohered distri- 
bution Q, but where the non-Gaussian observer will 
get, again, a highly squeezed Wigner distribution, with 
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-3-2-10123 



FIG. 7: The function w(u 2 ,0) normalized to its maximum 
value plotted against the physical coordinate normalized to 
the overall width of the distribution: f/y/NF, . The blue 
and red curves are the distribution inferred by the Gaussian 
and the non-Gaussian observers respectively for n = 10 (in 
the limit x 3> n 2 in the latter case). While the former sees 
a broad, i.e. incoherent, distribution, the latter instead finds 
highly squeezed distribution aroun ip 2 = ipg and, therefore, 
a highly coherent state. As already emphasized in Fig. 
the two peaks are in fact connected by the underlying O(N) 
symmetry and should not be interpreted as an incoherent su- 
perposition of two states. 



number of states whereas the extra piece of information 
the non-Gaussian observer has leads him to conclude in- 
stead to a state with a high degree of quantum coherence. 

An interesting extension of the present study will be to 
include other possible four-point correlators in the non- 
Gaussian description, such as (tt 4 ) or (n 2 ^ 2 ). In particu- 
lar, it will be interesting to see whether these change the 
present conclusions concerning entropy or purity. 

In the context of recent studies of entropy production 
and decoherence in quantum field theory based on the in- 
complete descritpion picture [a, |6| , the present work pro- 
vides a basis to study the role of dynamically generated 
non-Gaussian correlators concerning entropy production 
or the decoherence process. 

Finally, the present work is also of relevance for more 
general studies of nonequilibrium dynamics and thermal- 
ization in quantum field theory [12| . The existing litera- 
ture, see e.g. 0, li~7r{22j is, to its vast majority based on 
analyzing thermalization from a Gaussian observer per- 
spective, i.e. looking at the way equal-time two-point 
functions approach their equilibrium values [38j . We be- 
lieve our work opens the way for more detailed studies of 
the nonequilibrium flow in the space of correlation func- 
tions. 



5 V ~ 1 /n 2 <C 1 , thus exhibiting a high degree of quantum 
coherence. This situation can be observed for instance in 
the right-bottom (0 = it) plot of Fig. [6l which shows 
a slice in Af 2 -dimensional phase-space. Here the corre- 
sponding Gaussian distribution would be a long ellipse 
in the 7r-direction surrounding the two highly squeezed 
contour plots corresponding to the non-Gaussian distri- 
bution. 
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Appendix A: Connected four-point correlator C4 at 
large- N 



VI. CONCLUSION 

In this study we compare the effective density opera- 
tors inferred from two observers, one of which only mea- 
sures Gaussian correlators F, K and R and the other has 
further access to the field four-point function C4. As can 
be expected, we find that, depending on the strength of 
the measured nongaussiantity (quantified by our param- 
eter x), both observers can arrive at very similar or very 
different conclusions concerning various observables. 

Global observables such as the correlation entropy or 
the quantum purity are rather insensitive to the detailed 
structure of the underlying density matrix and thus to the 
amount of nongaussianity. In other words, both observers 
essentially agree about the degree of quantum purity of 
the system. 

In contrast, more detailed observables such as the 
probability distribution of field configurations, can be 
very different for sufficiently large nongaussianity. We 
identify situations where the Gaussian observer would 
conclude to a thermal-like, incoherent mixture of a large 



We derive the expression of the connected correlator C4 
in the limit N — > 00. The latter can be obtained in var- 
ious standard ways [T(| , for instance by computing In Z 
at next-to-leading order in the 1 /./V-expansion and using 
©, or by introducing a linearly coupled source term j a <Pa 
in the classical action and taking four functional deriva- 
tives of the resulting generating functional with respect 
to j, see e.g. [2j|. We present here an alternative, sim- 
ple derivation, based on introducing a bilinearly coupled 
source J(p a ip a - Consider the following functional (note 
that Z[J = 0] = Z) 

Z[J}= j VipV-Keti^^v^-Fivri+JvaVa}. (Al) 

J per. 

Using similar manipulations as before, it is an easy exer- 
cice to check that, up to an irrelevant constant, 



(A2) 
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Computing 51n Z[J]/5J(t)\j-q from either Eq. (|A1I) or 
Eq. (|A2|) . one gets the exact relation 



<X> = 2i V F, 



(A3) 



where 



pX F [x] ( TrLn G- 1 (x)+/ 1 dr X 2 /2r ) ) 

p xe -f(T r LnG-Mx)+/ 1 <ir X V2r ; ) ' (A4) 

Similarly, expressing dlnZ[0]/dB, one shows that 

F-2 J 4(TrG( X )). (A5) 

In particular, one has the exact relation 

( X ) =4iA V (TrG( X )), (A6) 

of which Eq. (TT4")) is the expression in the limit N —> oo. 
Next, from Eq. (|A"2j) . we get [11 



a " inz= -^ + ^r dT(x2(T)) - (at) 

Writing 

(X(r)x(r')) = -(2^) 2 + (x(r)x(r')) c , (A8) 

where (. . ,) c denotes the connected contribution and 
where we used (|A3|1 . and using Eq. (j6|), we obtain the 
exact relation 

N + 2 6(0) N 



2/'-' • — — f, 

N 2r) 



n r 1 

~ Wo dT(x2(T)) - (A9) 



Finally, we note that the connected x _ correlator can 
be obtained from 



1 8 2 \nZ[J] 



N SJ(t)6J(t>) 



= ^Ut-t>)-^- (x(r)x(r')) 
j=o 2f l V 2? 7 

(A10) 

We now consider the large- N limit. Up to an irrelevant 
additive constant, 



\nZ[J] 



--TrLnG-^x/)- ( * J + ^ , (All) 



where xj satisfies the saddle-point equation 
X.i +M = HArfTrGftj). 



(A12) 



Note that XJ=o = X: see Eq. (|14l) . Differentiating both 
sides with respect to J, one easily shows that 



8J(t') 



-iD(r — t'), 



(A13) 



j=o 



where the function £> is defined by (1 = 8(t — t')) 

D _1 = 1 + ??II with n=16A 2 G 2 (x). (A14) 



In particular, one has 

D = l-r)Il*D 



i-r7n + 77 2 n*n*L>, 



(A15) 



where [A*B](t,t') = f dr" A(t,t")B(t"t'). Differ- 
entiating Eq. (|A11I) twice with respect to J and using 
Eq. (|A10I) , one obtains the familiar expression of the x~ 
field correlator in the large- N limit fig : 



(A16) 



Inserting the latter in Eq. (|A9|) and using Trll = 
G 2 (x;0) = [TrG(x)] 2 = (2F) 2 , we get: 



G 4 



n 



Tr [n 7k- n * D] 



(A17) 



Appendix B: Matsubara sums 



Here, we compute the functional traces needed in Eqs. 
(Tl"4"| and (|A17|) , using standard methods of finite temper- 
ature field theory [2J] . The basic quantity is the Green's 
function (TT2"j) evaluated for a r-independent field config- 
uration x( T ) = X- Writing 

G(x;r-r') = .g(r-r',z), (Bl) 

where z = 2y/A(B' — ix), we have 

-g"(r,z) + z 2 g(T,z)=5{T). (B2) 

Periodic boundary conditions of the functional integral 
(0 results in the periodicity of the function g: g(r+l) = 
g(r) which can thus be written as a Fourier series: 



uj n = 2im. 



(B3) 



(B4) 



where 

Inverting (|B3|) . one obtains 

9(t, z) = ± ([1 + n(z)]e~' T + n(z)e^) , (B5) 
where n(z) — (e z — Eq. ([T5jl follows directly from 



TrG(x)=g(0,z) = J29n(z) = 



2ztanh(z/2)' 



(B6) 



The function IT = 16^4 2 G 2 (x), see Eq. (|A14|) . is easily 
obtained. Using the relations 



and 



i(z + z')[l + n(z) + n(z')] = n(z)n(z') (B7) 



n'(z) = -n{z)[l + n{z)] (B8) 
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one sees that the r-dependence of the square of the 
Green's function g(r, z) with "frequency" z is given by 
that of the same Green's function with frequency 2z: 



g 2 (T,z) = 2g(0,z)g(T,2z)- 



n'(z) 



(B9) 



Introducing 

6 n (z) = 32A 2 [ 5 (0, z) - n'(z)5 nfi ] (BIO) 

the Fourier components of the function IT read 

n»W = 9 n (z)g n (2z). (Bll) 

Finally, we write the Fourier coefficients of the \~ 
propagator D in Eq. (|A14|) as 



D- 1 (z) = l + V Il n (z). 



(B12) 



One easily checks that the r-dependence of D is essen- 
tially given by that of the Green's function g with a mod- 
ified "frequency" . Specifically, one has: 

U n (z)D n (z) = e n (z)g n (y{2zf + n6 n {z)) . (B13) 

The fact that, for what concerns their r-dependence, 
both II and IT * D are essentially given by the Green's 
function g allows us to rewrite the double-convolution 
in the expression (IA17P of C4 as a simple one-loop-like 
expression. For z solution of the saddle-point equation 
(TT5j). one has n(z) = n, i.e. z = ln(l + 1/n) and g(0, z) — 
1/2k. We obtain 



C 4 



>1 



, V g n (2z)g n (2z') + ? g»(2z)g {2z") 

(16kF 2 ) 2 [n^O 

(B14) 

where z' — z\J\ + x and z" = Z\J\ + (x. The sum is 
readily performed using [24[ 

1 r h( x ) + h(x') h(x) - h(x') 



^2 9n{x)g n {x') 



2xx' 



x + x' 



(B15) 

where h(x) = n(x) + 1/2. Equation (p~8|) follows, with 
yf(y) = h(2zy). 



with z the solution of the relevant saddle-point equation. 
We discuss the physically allowed values of z, i.e. those 
for which these integrals are well-defined. We consider 
only one O(N) component for the sake of the argument. 

The calculation of In Z involves an integral over fields 
with periodic boundary conditions ip(0) = <p(l) whereas 
(tp2\D\tpi) involves fixed boundary conditions <p(0) = tpi 
and <p(l) = if2- In the latter case, one can compute the 
explicit dependence in (pi and if2- 



Dip e 



-S[<p] = e -F u (z)u 2 -F„(z)v 2 



Dip e 



Slip] 



(C3) 



where u 2 = (ip 2 + </?i) 2 /4A, v 2 — (ip 2 — tpi) 2 /4:A and 
the functions F u (z) and F v (z) are given in Eqs. P5)l 
and Prj|) . The functional integral on the RHS is to be 
evaluated with Neuman boundary conditions: (p(Q) = 
<p(l) = 0. 

Writing ip(r) = J2 n e1,UJnT( Pn, one is finally led to eval- 
uate Gaussian functional integrals of the form 



Dtp e 



(C4) 



with different types of boundary conditions. Those are 
well-defined if w 2 , + z 2 > 0, Vn. 

For Neuman boundary conditions </?(0) = ip(l) = 0, 
one has u> n — mr with n G IN*. The lowest frequency 
in (IC4p is u>i = ir and the integral is well-defined for 

Z 2 > -7T 2 . 

For periodic boundary conditions <p(0) = <p(]~), the fre- 
quencies uj n = 2nn with n G Z and the lowest one is 
ujq = 0. The functional integral is well defined only for 
z 2 > 0. Alternatively, the integral with periodic bound- 
ary conditions can be obtained from the one with fixed 
boundary conditions as 



I D<pe~ s[v] = ( dtp [^Vpe 

J per. J J (p 



Slip] 



Dip e 



-Slip] 



dip e 



(C5) 



where we used (IC3[) . For real z 2 , the (^-integral is well- 
defined if F u (z) > 0, which is equivalent to z 2 > 0. 



Appendix C: On solutions of saddle-point equations 

The evaluation of both ln.Z and {ip2\D\ipi) in the limit 
N — > 00 involve Gaussian functional integrals of the form 



Dip e 



-Slip] 



(CI) 



with different types of boundary conditions, where 

<%] = A fdrip (C2) 



4A 



Appendix D: Wigner function and quantum 
coherence 

Here we briefly recall the relation between the shape of 
the Wigner distribution on phase-space and the degree of 
quantum coherence between distinct elementary cells of 
phase-space, described by coherent states. The latter are 
the eigenstates of the usual annihilation operator [4pj : 
a a \a) = a a \a) ,Va, with 



V2 



(Dl) 
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Recall that the Wigner distribution corresponding 
to a given coherent state is a Gaussian centered in 
\/2(Rea a , Ima a ) and of width l/v2 in both ip and tt 
directions, in units of the natural field dimension (3. [25|. 
Because in that state the quantum fluctuations in both ip 
and 7r are as small as they can be according to the Heisen- 
berg uncertainty principle, it is often called semi-classical 
states. It can be seen as describing an elementary quan- 
tum cell in phase-space. 

Using 



V2 



and 



V2 



-+q -a 



(D2) 



(D3) 



where \a\ 2 = a* • a, one easily derives the following for- 
mula for the matrix elements of a given density matrix 
D in terms of the corresponding Wigner function: 

(a\oc) J 7r JV 

(D4) 

where p v = and & = . 

Let us first recall the simple case of a Gaussian Wigner 
function centered at the origin |4lj : 



in the case of interest here, see Eq. (|71l) with F/n 2 <C 1, 
is given by (p* ipo- Depending on the value of in 



(|D8|) . one may employ different approximations for the 
Bessel function and evaluate the integral. For instance, 
for |/3 v <po | 3> N/2 one can use the leading asymptotic 
behavior J u (z) ~ \J2~fiTZ cos(z — vtt/2 — 7r/4) and is then 
left with the evaluation of Gaussian integrals. We obtain 
finally, in terms of real and imaginary parts of a and a' : 



\{a\D\a') 



l^y) |co S (2i/3^o-AT7r/4)| 



(g+gOj (S+S') 2 

2(l + 2Aj) 2(l+2Aj) 



')2 A 2 (b-?) 2 



x e 



(D9) 



As an illustration, consider two phase-space cells such 
that d — a'* — (ip e + ix)/y/2. One gets: 



\(a\D\a')\<x 



pl + x 2 



A cp 2 
-T-^-X 

'A- 



(D10) 



The conclusions drawn from the previous case, Eq. (|D6I) . 
concerning the degree of quantum coherence still hold. 



Appendix E: Peak structure of d{u 2 



W(tp, tt) (X e 2& f 



(D5) 



for which the calculation is straightforward. In terms of 
real and imaginary parts a = a + ib and a' = a! + ib', 
one obtains 



We discuss the structure of the matrix element (|5Tj) for 
n and x sufficiently large so that a clear peak structure 
appears, see e.g. Fig. [3] For n 3> 1 and a; 3> 1 the latter 
is well-described by a Gaussian. Recalling Eqs. ([57)1 and 
we have 



(a\DW) 



(s+a') 2 



(S+S') 2 



\ 2 i^ — Fj') 2 



^|(S-S') 2 



xe 



(D6) 
quantum 



One thus observes non-trivial correlations - 
coherence - between distant phase-space cells 
1) when the Wigner distribution is sufficiently squeezed 
in either the tt (A t < 1) or the tp (A y < 1) direction. 
Note that the size of the correlation in the ip (resp. tt) 
direction - i.e. along the real (resp. imaginary) axis in 
the coherent state basis - is solely controlled by the width 
of the Wigner distribution in the 7r (resp. ip) direction. 

Consider now the case of interest in Eq. (|69| : 



W{0, tt) oc e 



(D7) 



The 7r-integration goes as in the previous case. As for 
the (^-integration, N — 1 angular integrals are easily per- 
formed and one is left with the following radial integral: 



dtpip^e 2A ii JN^i^cp), (D8) 



d u hid 
d v \nd 



-2uF u (z) 
-2vF v (z), 



(El) 
(E2) 



where z = z(u 2 



Since F v (z) > for real z 2 , the only 



extremum of hid in the u-direction is located at v = 0. 
It is a maximum. Instead, d u Ind vanishes both at it = 
and at the point where F u (z) = z(u 2 ,v 2 ) = 0, i.e. 
u 2 = u 2 (v 2 ) see Eq. (|59|) . The latter extremum only 
exists if conditions (|53")l are satisfied, in which case it is a 
maximum whereas the one at u — is a local minimum. 
We consider this case in the following. 

The absolute maximum of the function lnd in the 
(u, u)-plane is therefore located at (uq,0) with 



1 



(E3) 



Performing a Taylor expansion near this point (denoted 
by a subscript 0, one has 



d 2 u \nd\ = -4u 2 F^z)d u 2z\ 



(E4) 



with P v = V ■ /3 y ) 1/2 . For N > 1 the factor before the 
Bessel function is strongly peaked around ip = tp* which, 



d 2 lnd\ a = -2F V (0) 



(E5) 
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and the cross-derivative d u d v In d| = 0. Prom the saddle- 
point equation (j4"T|). one gets 



d„2z — 



fu(z) 



(E6) 



Using small- z behaviors of the functions Fi(z) and fi(z), 
with i = 0, u, v, we get, after calculations, 



In d = In do 



(u - u ) 2 _ %P_ 
25 2 251 



(E7) 



where do = d(u^,Q), with 



{ «=4G + ir + i) md 45 4 (E8) 

For ra> 1, one has Mg ~ 2n 2 (l — l/2x), £ ?z x/n 4 and 
O (n- 2 ) and ^ = i (E9) 



2x- 1 6 



The result of main text follows in the limit x ^> n 2 . 

Finally, we note that a similar calculation yields, for 
the first non-vanishing higher-order derivatives, in the 
same limit. 



d 3 u lnd\ 



3u 5l \5l 



1 



and 



d 4 \nd\ a = 



3u 2 5 2 u 



(E10) 



^^i^GH^ " 1 ,E| " 



(E12) 



The function d is thus well-described by a Gaussian, 
which is confirmed by our numerical results. 
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